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We consider the problem of the renormalization of the vacuum polarization in a symmetry space- 
time with axial but not spherical symmetry, Schwarzschild space-time threaded by an infinite straight 
cosmic string. Unlike previous calculations, our framework to compute the renormalized vacuum 
polarization does not rely on special properties of Legendre functions, but rather has been developed 
in a way that we expect to be applicable to Kerr space-time. 



I. INTRODUCTION 

o 

f^ In this paper we calculate the vacuum polarization of a massless, minimally coupled scalar field in the region exterior 

OQ to the horizon of a Schwarzschild black hole threaded by an infinite thin cosmic string which reduces the spherical 

(-H symmetry to axial symmetry. We consider the scalar field in the Hartle-Hawking vacuum state, corresponding to a 

^ black hole of mass M in (unstable) thermal equilibrium with a bath of blackbody radiation. 

^~~i There is an extensive body of work on renormalization of the vacuum polarization on black hole spacetimes ( [1- 

C^ 5]). In all of these cases, the authors have considered spherically symmetric black holes. The most astrophysically 

C^ significant case, however, is the Kerr-Newman black hole. The calculation of {(p'^)ren in this case has proved elusive 

(with the exception of it's calculation on the pole of the horizon where the effects of rotation are minimized [6]). The 

fj principal reason for this is that the rotation of the black hole means that it no longer possesses spherical symmetry, 

^^ but only axial symmetry. The existing calculations, referenced above, rely heavily on the Legendre Addition Theorem 

' and other well known properties of the special functions, as well as the Watson- Sommerfeld formula. In the axial 

(^ symmetric case, the Legendre Addition Theorem does not apply and the Watson-Sommerfeld formula is no longer 

'— ^ useful as it is no longer possible to perform the sum over the azimuthal quantum number. 

, In the case being considered in this paper, we are dealing with a black hole where the symmetry is reduced from 
spherical to axial but without the added complication of rotation. This fact makes this case an ideal precursor to the 
Kerr-Newman case. It presents the first calculation of the renormalized vacuum polarization on the exterior region 

l/~) of an axially symmetric black hole. It should be noted that DeBenedictis [7] has calculated the vacuum polarization 

^^ on an axially symmetric metric, but in the case of a black string, not a black hole. Most importantly, the method 

^^ presented in that paper is not applicable to the current space-time or to the Kerr-Newman case. 

r — The method we present here does not rely on specific properties of the angular functions (such as Addition Theo- 

(^ rems) and that results in a complete modc-by-mode subtraction (as opposed to the partial mode-subtractions in the 

^^ literature). Furthermore, we elucidate some important points about the Christcnsen-DeWitt point-splitting approach 

^~~^ [8] to renormalization. In particular, we show that the choice of point-separation direction is intimately connected to 

^ the order of summation of the mode-sum. 
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II. THE MODE-SUM EXPRESSION FOR THE GREEN'S FUNCTION 

The cosmic string is modeled by introducing an azimuthal deficit parameter, a, into the standard Schwarzschild 
metric. We may describe this in coordinates (t,r,0,(l)), where 4> is periodic with period 2TTa, so we may take (j) G 
[0, 2Tra), in which the line element is given by 

ds^ = -(1 - 2M/r)dt^ + (1 - 2M/r)"Mr2 

+r^d9^ +r^sm'^9d^^. (2.1) 

For GUT scale cosmic strings a = 1 — 4/x where /z is the mass per unit length of the string so we shall assume 
< a < 1. In Sec. V, we further limit ourselves to the case 1/2 < a < 1 which is physically justifiable since /i ^ 1. 
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We can alternatively define a new azimuthal coordinate by 

^ = a(t) (2.2) 

so that (j) is periodic with period 27r and we may take (p G [0, 27r). In coordinates (i, r, 9, 0), the line clement is given 

by 

ds^ = -(1 - 2M/r)dt^ + (1 - 2M/r)-^dr^ 

+r'^Ae'^ + a^r'^ii\T?eA4P. (2.3) 

We shall consider a massless, minimally coupled scalar field, ip, in the Hartle-Hawking vacuum state. Since this 
is a thermal state, it is convenient to work with the Euclidean Green's function, performing a Wick rotation of the 
temporal coordinate t — > —it and eliminating the conical singularity at r = 2M by making r periodic with period 
2tt/ K where n = 1/(4M) is the surface gravity of the black hole. The massless, minimally coupled scalar field, satisfies 
the homogenous wave-equation 

n^(T,r,0,0) = O, (2.4) 

which can be solved by a separation of variables by writing 

ip{T, r, e, 0) - e"''^+^"'^P(^)i?(r) (2.5) 

where P{0) is regular and satisfies 

[9] = (2.6) 
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The A(A + 1) term arises as the separation constant. The choice of A is arbitrary for ip to satisfy the wave equation 
but requires a specific choice in order for the mode-function P{9) to satisfy the boundary conditions of regularity 
on the poles. In the Schwarzschild case (no cosmic string, a = 1), regularity on the poles means that X = I, i.e the 
separation constant is l{l + 1). In the cosmic string case, the appropriate choice of A that guarantees regularity of the 
angular functions on the poles is 

X = l-\m\ + \m\/a. (2.8) 

With this choice of A, the angular function is the Legendre function of both non-integer order and non-integer 
degree, wiz., 

P(6i) = P^"l"l''"(cos6i). (2.9) 

It can be shown that these angular functions satisfy the following normalization condition, 

/rJi':'fi™i/.(cos^)p^7j/;„,/jcose)d(cos^) 
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The periodicity of the Green's function with respect to {t~t') and [cfy — cj)') with periodicity 27r/K and 27r, respectively, 
combined with Eq.(2.10) allow us to write the mode-sum expression for the Green's function as 



G(-'-') = £ E ^"''^^"^'^ E e-(^-^')^(2A + l) ^|^+[;;;[/;^ + |j p-l"'l/"(cosg)P,--l"l/°(cosgOXnA(r,/). 

(2.11) 
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where Xn\{r,r') satisfies the inhomogeneous equation, 

{d;:(^^-2M.)^-A(A + l)-^:^-^}x„.(r,r') 

= --(5(r-r'). (2.12) 

a 

It is convenient to write the radial equation in terms of a new radial variable r\ — r jM — 1, the radial equation then 
reads 

= -\m^-tI). (2.13) 

where we have used the fact that k — 1/(4M). For n = 0, the two solutions of the homogeneous equation are the 
Legendre functions of the first and second kind. For n 7^ 0, the homogeneous equation cannot be solved in terms of 
known functions and must be solved numerically. We denote the two solutions that are regular on the horizon and 
infinity (or some outer boundary) by Pnxi'']) and q-nxi'']), respectively. A near-horizon Frobenius analysis for n 7^ 
shows that the indicial exponent is ±|n|/2, and so we have the following asymptotic forms: 

(2.14) 

Defining the normalizations by these asymptotic forms and using the Wronskian conditions one can obtain the 
appropriate normalization of the Green's function: 

, \fPx{v<)Qxiv>) n-0 
X«A(^,r/')-<^"^l (2.15) 



Pnxiv) " 


^(77-l)l"l/2 


V- 


-> 1 


Qnxiv) " 


. (^_l)-|n|/2 


77- 


> 1. 



01 I njPnx{n<)Qnxiv>) n^O. 
2\n\aM 

III. CHOOSING A SEPARATION DIRECTION 

In order to renormalize (<^^), we subtract, in a meaningful way, the geometrical Christensen-DeWitt renormalization 
terms from the Green's function. This approach rests on the fact that all Green's functions will possess the same 
short distance behaviour encapsulated in the Hadamard parametrix. On the other hand the full Green's function 
must refiect the relevant boundary conditions of the global problem, for example periodicity with particular period in 
T and (j), and this is most easily expressed by a mode decomposition. In order to perform the renormalization it is first 
necessary to perform a regularization of the Green's function and the natural rcgularization in the Christensen-DeWitt 
approach is to consider the Green's function with the two points separated. The geometrical singularity may then be 
removed prior to bringing the two points together to give {(fP').ren- The geometrical nature of the subtraction means 
that we obtain the same result whichever direction we separate in, although there is a surprising twist in the tale 
described below. 

As we are free to choose the direction of separation, we may do so in a way which makes the calculation as 
straightforward as possible. In almost all black hole calculations in the literature, e.g. [2-5, 9], the authors have 
preferred to separate in the temporal direction, except for on-horizon calculations. The principal reason for this is 
that the metric components do not depend on t, making the renormalization terms somewhat easier. (The same 
could be said for separating in 0, however, making it an equally suitable candidate.) For on-horizon calculations, 
it is typically most convenient to separate in the radial direction (see [1] for example). In fact, we have calculated 
elsewhere [10] an analytic expression for {(fP')ren on the horizon of the Schwarzschild black hole threaded by an infinite 
cosmic string by separating radially and using a summation formula we have derived in that paper. 

Since we have calculated the vacuum polarization on the horizon in another paper, we shall concentrate here on 
the calculation of off-horizon values in this paper. The question we must address is what is the most convenient 
separation when we no longer have spherical symmetry. To answer this we start by analysing the approach taken in 
the literature in the spherically symmetric case. In this case, one can sum over m to obtain expression for {(p^)ren 
that involves an inner sum over /-modes and an outer sum over n-modes. One then converts the /-sum into an integral 
using the Watson-Sommerfeld formula and sums the n-modes directly. (This is, of course, only a sketch of the method 



and several other tricks and techniques are used to do the calculation, none of which are important to the choice of 
separation.) 

In the cosmic string case, the axial symmetry means that converting an Z-sum to an integral results in numerical 
integrals over the square of Legendre functions; this is neither convenient nor useful. The most practical way to 
proceed with the calculation is to convert the n-sum to an integral and sum the I, ?7i-niodcs directly, since the angular 
part of the Green's function does not depend on the n-modcs. In fact, this proves to be a very fruitful approach to 
these calculations, both in the axially symmetric and the spherically symmetric case. There is a caveat, however, 
which relates to the order in which we perform the mode sums which is intimately related to the distributional nature 
of the expressions we are dealing with. 

At a 4-dimensional level we have a sum over a complete set of mode functions and the expression 
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is understood in the sense of smearing with smooth functions of compact support in the 4-dimensional space. In 
moving to a point separated expression with just one coordinate different we must consider the S convergent limits in 
three of our coordinate directions. In particular as d/dr and d/d(j) are commuting Killing vectors we may associate 
with them independent quantum numbers n and m. Also correspondingly, when we derive our Green's function, it 
matters not whether we separate out the t dependence or the dependence first. However, when we limit our test 
functions to 3-dimensional delta functions the corresponding sum must remain until last, i.e., the outer sum must be 
that which corresponds to the direction in which we have separated. More specifically, separating in the temporal 
direction corresponds to an inner Z, m-sum and an outer n-sum, separating in the azimuthal direction corresponds to 
an inner n-sum and an outer /, m-sum. 

We have two strong argumets in support of the above argumet. Firstly, one can show analytically, in the 
Schwarzschild case, that the finite difference between summing in the different orders is precisely the analytic differ- 
ence between the regular parts of the Christensen subtraction terms for temporal and azimuthal separation. In other 
words, separating in the temporal direction and doing the n-sum first gives a finite but wrong answer! An alterna- 
tive way of understand this is that the renormalization procedure adopted by Candelas-Howard and the alternative 
procedure of this paper results in a double mode-sum that is convergent, but not absolutely convergent. As a result 
the order of summation matters. 

The second supporting argument is that we have an unambiguous answer on the horizon in both the Schwarzschild 
black hole with a cosmic string [10] and without a string [1]. Clearly the correct off- horizon answer must match up 
with the horizon value as the horizon is approached and this is only true when we sum according to the rules laid out 
above. 

It is natural for us to separate in the azimuthal direction for the cosmic string case and so in light of the previous 
arguments, the appropriate mode-sum form for our unrenormalized Green's function is 
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(3.1) 
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IV. RENORMALIZATION 



For a massless scalar field in a Ricci-fiat space-time, the only Christensen-DeWitt subtraction term required for the 



calculation of {ip'^ 



IS 



^div \'^j ^ ) 
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where 2a{x,x') is the square of the geodesic distance between x and x' [8]. For an azimuthal splitting this becomes 
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One is now faced with the challenge of subtracting this geometrical expression from the mode sum Eq.(3.1) in such a 
way that the limit may be performed. The approach in the literature, e.g. [2-5], is to bring the divergent term inside 
the outer sum (in this case the m-sum, in the temporal splitting case the n-sum) using an identity from distribution 
theory, 
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or its equivalent expression for temporal separation. The success of this approach in the spherically symmetric case 
relies heavily on the applicability of the Legendre Addition Theorem and the Watson-Sommerfcld formula. This is 
not possible in the general case. Instead we would like to be able to write the divergent term in Gdiv as a triple 
mode-sum over m, I, n so that a full mode-by-mode subtraction may be performed. 

It turns out that there is a natural and general way to approach this problem. There are a whole gamut of summation 
formulae that can be derived simply by equating different but equivalent expressions for the same Green's function. 
This requires little knowledge of the special functions being summed, only that they are solutions to the homogeneous 
wave equation being considered. Indeed, this approach is one of the standard ways of proving the Legendre Addition 
Theorem. The universality of the geometrical singularity structure provided by the Hadamard form then ensures, 
that with appropriate parameterisation, we can match the required coordinate divergence to an appropriate mode 
sum. 

Ideally, we equate a mode-sum expression for a Green's function to a closed-form or quasi-closed form expression. 
The most effective way of deriving such summation formulae is by considering appropriate Green's functions on 
Minkowski spacetime, where the Green's function is usually known in closed form or quasi-closed form. In the next 
section, we will derive the appropriate summation formula for the current problem by considering the thermal Green's 
function on Minkowski spacetime threaded by an infinite cosmic string. 

V. THERMAL GREEN'S FUNCTION ON MINKOWSKI SPACETIME WITH A COSMIC STRING 

We consider a massless scalar field at non-zero temperature T propagating in Minkowski spacetime with an infinite 
cosmic string running along the polar axis. The Euclideanized metric is given by 

ds^ = dr^ + df + ^^d9^ + a^f sin^ 6d(j)'^ , (5.1) 

where 9 and (p are the usual polar coordinates on the 2-sphere and £, is an arbitrary radial variable. The periodic 
(thermal) Euclidean Green's function G^(At, Ax) can be written as an image sum over the zero temperature Green's 
function, G'(At, Ax') as 

oo 

G0(At,Ax)= Y^ G'(At + fc/3. Ax) (5.2) 

fc— — oo 

where j3 = 1/T is the inverse temperature. A mode-sum expression for the scalar field Green's function at zero 
temperature is easily found to be 

G(Ar,Ax) . -V Te-^Mc. f e^^^ f (2A + 1) ^}^ + H^" + jj 

p-|"l/"(cOS0)p-|"l/"(cOS0')(^^/A+l/2(^e<)i^A+l/2(c^C>) (5.3) 

where / and K are the modified Bessel functions [11] of the first and second kind, respectively. Now, using the Fourier 
transform of a Comb function, 

E e^-'^^'^-j £ 5{u-n2n/Pl (5.4) 
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we arrive at an appropriate expression for the thermal Green's function: 
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where n — 2ttT, and the n = term is understood in the sense 

hm(ee')"'/'^A+i/2(ri<<)i^A+i/2(ri<>) = (^)' g^(2A + l) - ^^'^^ 

In particular, for azimuthal separation, we have (taking into account that we must choose the appropriate order of 
summation) 

m=-oo l=\m\ V I 1/ ; 

oo _, 
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Note that, with the exception of the radial functions, this is completely equivalent to the Green's function Eq (3.1). 
We have an equivalent expression for this Green's function given by Linct [12], valid for 1/2 < a < 1, in which 
the author writes the Green's function as the singular part plus a regular integral part, G^(x,a;') — Ggingi^T^') + 
Gi„t{x,x'), The integral part Gmt is given by 

G,nt{x, x') = -/(Ac/. + 7r/a) + /(A0 - n/a) (5.8) 

where 

m) = ^^^^(^) r siuHnRju)) 1 ^^ 

■'^ ' Sn^a Jo R{u)[cosh{KR{u))-cosKAT]{cosh{u/a)~cos{^)) ^'^ 

and 

i?(u) = [C^+C'^-2^^'cos0cos0' + 2CC'sin6lsin6i'coshu]^/^ (5.10) 

The singular part is 

Gsing\X,X ) ~ —^ (5.11) 

4tt p[cosh up — COS kAt\ 

where p = [C^ + C'^ - 2^^' cos 9 cos 9' - 2^^ sin 9 sin 9' cos(a A</-)] ^/^ . 
In particular, for azimuthal point separation, we have 



G^(e, 9, A0) = ,_,,,^,^,„2^^^, + T^^T^r^^ + T^ + G^nt{^, e, Acf,) + O(A0)2. (5.12) 
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Finally, equating Eqs. (5.7) and (5.12), we arrive at the very useful identity 
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It is important to emphasize that this equation is true for any T and for all ^. For our purposes, the obvious choice for 
T is the temperature of the Schwarzschild black hole so that k is now the Schwarzschild surface gravity k ~ 1/(4M). 
In the next section we show that there is a prescription to assign ^ in terms of the Schwarzschild radial variable rj in 
such a way as to guarantee the convergence of the mode-sum. 



VI. MODE BY MODE SUBTRACTION 

We now return to the Schwarzschild cosmic string case. From Eq. (3.1) and Eq. (4.2), we have the following 
expression for the renormalized vacuum polarization: 

(^2)^^„ = lim [G(r, 9, A0) - Gd^v{r, 9, A,/.)] 
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Now our identity, Eq.(5.13), is of precisely the correct form to allow us to convert the 1/A02 term into an appropriate 
triple mode-sum. On dividing Eq. (5.13) by M^{r] + 1)2 and substituting into Eq. (6.1), we obtain 
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We can now take the limit inside the sum to get 

(<^ )ren = (<^ )s«m + (<^ ) int + (<P )analytic, (6.3) 

where 
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where i?(u) = (2^2 §1^2 ^^^ + coshu))^/2_ •\Ye have used Eqs. (5.8) and (5.9) to obtain the expression for {(p^)int- We 
have also explicitly separated out the n — term. 

A key feature of our approach is that since we have a triple mode-by-mode subtraction, with the correct choice of ^, 
we can ensure the convergence of the triple sum, thus making redundant the removing of the 'superficial' divergence 
discussed by some authors ([3, 5]). Specifically, we choose ^, such that, for a given n, the I, m mode-sum is regular. 
We can do this by associating the /, m mode-sum with a particular 3D Green's function by a technique known as 
dimensional reduction [5], and then we examine the Hadamard singularity structure of this Green's function. In the 
Appendix, we show that the condition of regularity of the I, m sum can be achieved by taking ^ to be 

^ = ^^ (^2_i)i/2 - (6-7) 

The mode-sum now converges for this choice of ^, as we discuss in detail in the following section. 



VII. WKB APPROXIMATIONS 

In order to analyze the convergence of the mode-sum of Eq. (6.4), we consider the WKB approximations of both 
the radial part of the Green's function and the subtraction terms. We have adapted the WKB prescription given by 
Howard [2, 5], valid for large n and I. Employing a second order WKB approximation, we have 
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where /3„;j , f^l^x ^^'^ PnX ^^^ ^^^ zeroth, first and second order approximants respectively. Each term is written in 
successive powers of 'in\{'n)~^i where 
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and ujn — {n/4){'q +1)^. In terms of ^„Ai the approximants are 
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We also require the WKB approximation to the subtraction 
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Immediately, we see that the zeroth order terms are equal and therefore the slowest order term in the mode-sum is 
proportional to ^~^. From Eq.(7.2), this implies that for large I and n, the summand is 0{1/{P + n^)'^/^), since the 
angular functions scale linearly with /. Thus, it is now clear from the cancellation of the zeroth order approximations 
that the mode-sum converges. This proof of convergence is considerably simpler than analagous proofs in the standard 
approach (e.g., [2, 5] ). However, though we have shown that the mode-sum converges, the convergence is extremely 
slow (as is usually the case with such calculations) and is not absolute as discussed in Sec. HI. 

A standard trick for speeding the convergence in order to make the mode-sum calculation amenable is to subtract 
and add the second order WKB approximations [2]. For the n = terms, the subtraction term is exactly the zeroth 

order approximant to the radial Green's function and so we need only subtract and add (3^^ and /3q_^ terms in this 
case. We now have the following expression for the mode-sum 
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where A is the slowly converging part of the mode-sum. We further write 

S = El - S2 + I/2E3 
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n^l /^O 771— — I 

oo / 



:^ _V V roA I 1^ ^(A + |m|/Q + l) |„l/„ 2f PA(r?)QA(?7) .(o) .(i) .(2)1 



i=0 m=-i 



(7.^ 



Wc have re- written the I, m sum and swapped the order of summation with the n-sum which presents no problem for 
these sums since they are aU rapidly and absolutely convergent. In fact, the summand of Si and S2 are 0{l/{P+n^y^^) 
for large I and n. The summand of E3 is 0{l~^) for large I. 
The slowly convergent term, A, is given by 



(1) fl(2) 



* - s E (^A + ') m-!:::!/:t!! ^^''""'°'°"''>'{¥ + ¥ + s (ft- - -2 + a- - -s) } 



1=0 m=-l 



(7.9) 



In order to speed the convergence of this sum, we can convert the n-sum here to an integral using a modified version 
of the Plana- Abel Sum Formula. 



VIII. MODIFIED PLANA-ABEL SUM FORMULA 

We begin by re- writing the difference of the WKB approximants of order i (i = 1, 2) as [13] 



Pn\ InX Z^ („+ l)4»+2 r^2 _^ „21»+j + l/2 

j— ^ ' ' L J 



where f2 is given by 



n 



4(A + l/2)(r72- 1)1/2 



iv+iy 



and the coeflficients Cij{ri) are tabulated in Table 1. 



(8.1) 



(8.2) 



TABLE I: The coefficients Cij(r]) of the difference of the WKB approximants. 



c^An) 





1 


2 


3 


4 


1 


1 

16M 


81? (2') -3) 


I6m(4»7-5) 






2 


25Im(16'?' + 11) 


^(12,73-24,72 
+ 14,7-31) 


1^(6^-177,7' 
-^204,7 - 122) 


i^(lV-43,2 
+52,7 - 26) 


2if(8^'-26^' 
-h32,7 - 15) 



We now convert the ri-sum above to an integral using a modification of the Plana- Abel sum formula [13, 14]. In 
our particular case, we have 



E 



,2^ 



5n\ 



JO 



^^ [ii^ + n^Y+i+y^ 2r22i+i 7o [n^ + n^y+j+i/ 



,2J 



:— dn + 2(-l)*+^- 



where 



H-s) = 



(5e«^/2)2j 



(S -I- f])i+J + l/2 £2^^^ - 1 ■ 



r(. + j + 1/2)7^ (.-17)1/2^^ (8.3) 



(8.4) 
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Furthermore, the n-integration of Eq.(8.3) can be done exphcitly using 



"''■ dn=^^^:)^((+^/^). (8.5) 



/o [f72+n2]'+J+i/2 n^'2T{i + j + l/2) 

Substituting Eq.(8.5) into Eq.(8.3), we get 

V..W_ (.)^ 42'+iQo(r;) 42^+1 ^ Q,(r;) r(»)r(j + 1/2) 

Z^ IP„A ln\) 2(77 + l)4»+202»+l (?7 + l)4»+2 ^ 172^ 2r(z + J + 1/2) 

H — i J — u 

42'+i ^ 2a,iv)i-iy+^v^ /- /.('+■'•)(.) 

+ (^+l)4.+2Z. r(7+J + l/2) ^ (S- 0)1/2 ^-^ 

Using the coefficients of Table 1, the first term here is conveniently just — /3oa/2, while the second term here vanishes 
for i = 1,2, that is, 

-a(,)r(, + i/2) ^^^ ^^ ^^ 

^t^ r(z + j- + i/2) ^ ^ 

Finally, we arrive at the following expression for the n-sum of the difference of the WKB approximants of i*'' order: 

V (fiW - .«^ - -^ + 42-+1 ^ 2Q,-(r;)(-l)-+^-7^ r°° Z^'^+^'H^) , .„ „. 

Z. l/^nA 7„AJ - 2 ^ (r; + 1)4^+2 ^ r(z + J + 1/2) io (s - 0)1/2 ^* l»-»^ 

n — i J — U 

Returning now to our expression for A, we substitute (8.8) into (7.9). Remarkably, the n = terms cancel, leaving 
only the rapidly convergent integrals of Eq.(8.8), viz.. 



2-^.i-J ^n>^-\-^\/-+^) ' ^ ^^hh(v+ir^^' r(z+, + i/2) L j^^w^^^'r 

(8.9) 



m= — 00 l = \m\ i=l j=0 



IX. NUMERICAL EVALUATION OF {<^'^) 



c2\ 

ren 



We have shown that we can write the renormalized vacuum polarization outside the horizon as an analytic part 
Eq. (6.6), a regular integral part Eq. (6.5) and a contribution coming from the mode-sum (Eqs. (7.6)-(7.8) and 
Eq. (8.9)). In addition, we have calculated elsewhere [10] the vacuum polarization as an analytic expression on the 
horizon. 



,'i2\horizon / 1 i V 



Combining these results gives the renormalized vacuum polarization on the entire spacetime region of interest. 

We now turn to a specific example. Thus far, we have described the calculation of {(p'^)ren on the exterior of 
the Schwarzschild spacetime threaded by a cosmic string from the horizon out to infinity. The results we shall 
show, however, are for a Schwarzschild blackhole threaded by an infinite cosmic string inside a spherical box. The 
reasons for this are twofold: Firstly, it is numerically easier. Rather than evaluating numerical radial modes over 
the entire radial range, we are now limited to a finite range between the horizon and the boundary. Secondly, this 
calculation is motivated by the analagous Kerr renormalization. In this case, one cannot define a Hartle-Hawking 
vacuum everywhere on the exterior of the black hole. However, in order to determine the vacuum polarization on a 
state that possesses the defining features of a Hartle-Hawking vacuum, one must put in a spherical mirror inside the 
speed of light surface [15]. In any case, as a check of our method, we have also run the calculation where we take the 
boundary out to a very large radius, which yields the Candelas-Howard result [2] in the a — > 1 limit. 

For the results shown, we have taken units where the mass M ~ 1, the azimuthal deficit a — 0.95 and the boundary 
imposes Dirichlet boundary conditions at r^f, = 3 (or r^ = AM). We take an (r, 6) grid that consists of 50 radial points 



11 




<0'>,e 



FIG. 1: Plot of (87rAf)^{(^^),.en- It is clear from the figure that the vacuum polarization diverges at the poles, as we would 
expect since there is a curvature singularity there. 




FIG. 2: Plot of {STvM)^{ip^)ren for particular values of cos^. We observe that the effect of the string only becomes important 
very close to the string. This is due to a being close to 1. 

and 70 angular points. This boundary condition requires that the Green's function vanish on the mirror, implying that 
the outer radial function becomes 



P7iX\'lb) 



(9.2) 



On the horizon, we now have 



if )r, 



orizo7i,box 
ren 



1 



1927r2Af2 



cP- sin" B 



1 



OO / 



327r2Af2a ^-^ ^-^ 

1=0 m=-l 



(2A + 1) 



r(A 



r(A 



\m\/a ■ 



-4p-|'"l/"(cos0)2' 

1) ^ ^ ' 



<iVb) 



PxiVb) ■ 
(9.3) 
The sum here converges rapidly (about 12 decimal places after only 20 ^modes). 

With the exception of the boundary terms above, the calculation inside a spherical box remains unchanged. The 
calculation of {'p'^) analytic is trivial and requires no discussion. Similarly, the calculation of {ifP')int is numerically 
evaluated over the grid without difficulty. As we have mentioned the sums of S (Eq.(7.8)) are 0{l/{n^ + Py^^) and 
are therefore rapidly convergent. Close to the horizon and the boundary, however, this behaviour doesn't become 
apparent in Ei and E3 until the higher I and n-modes, compared with the interior region of interest. For this 
reason, after 100 Z-modes and 20 n-modes, we see convergence to about 6 decimal places near the horizon and the 
boundary and about 10 decimal places in the interior for Ei and E2. We see convergence to at least 10 decimal places 
everywhere for E2 for 100 Z-modes and 20 n-modes. The expression A (Eq.(8.9)) is exponentially convergent. Again, 
this convergence does not become apparent until the higher Z-modes near the horizon. We observe convergence to at 
least 10 decimal places everywhere in the region of interest by taking 30 Z-modes near the horizon and 10 Z-modes 
away from the horizon. 

We see from the graphs that the Casimir effect, arising from the presence of the boundary at rji, = 3, actually 
dominates when the boundary is this close to the horizon and becomes divergent as we approach the boundary. We 
also note that the vacuum polarization increases (becomes more positive) as we get closer to the poles. Again, we 
would expect this since there is a curvature singularity on the poles due to the cosmic string. In this particular case, 
since the azimuthal deficit is so small, the effect of the cosmic string is almost negligible until we get very close to 
the poles, responsible for the sharp divergence in Fig. 1. This graph would be more rounded and smooth for smaller 
values of a (corresponding to larger azimuthal deficits). We see from Fig. 2 that the cosmic string has little effect on 



12 




FIG. 3: Plot of (SvrAf )^((^^)ren for particular values of cos( 




FIG. 4: Plot of (SvrAf )^((^^)ren for particular values of cos( 



the vacuum polarization when we are not close to the poles. The 9 dependence diminishes the further we are from 
the horizon, i.e., when we are not close to the pole of the black hole, the effect of the cosmic string diminishes as we 
move further from the horizon. This is evident by comparing Fig. 3 and Fig. 4, since we can see the merging of the 
graphs as we move away from the horizon. The Casimir divergence due to the boundary is also clear from Fig. 4. 

The dominance of the Casimir effect serves to obscure the behaviour one would expect to see in the absence of a 
boundary. However, we may subtract the Casimir divergence by doing a DeWitt expansion of the heat kernel about 
the boundary [16]. For our particular geometry of spacetime and surafce, we find that the expected divergence from 
the curved boundary is given by 



m 



2\C asimir 



' ( ' 



1 



327r2M2 V {rib - rif I2{r]b - v) 



(9.4) 



Subtracting these divergent terms provides a clear picture (modulo finite boundary effects) of {'f>^)ren hi^ the absence 
of a boundary without evaluating the mode-sum over the entire exterior region of the black hole. In particular. Fig. 5 
possesses the general features of the renormalized vacuum polarization without a boundary. This is very evident 
from Fig. 6 where we have plotted {(p'^)ren for the Schwarzschild case (a — > 1, i.e. no cosmic string) using the data 
tabulated by Candelas and Howard [2] with the plot of {ip^)ren — (<^^)di"*'*™*'^ ^^^ ^^^ equatorial plane. The profiles 
are clearly identical, differing only by a constant finite amount, though this finite difference is constant only in the 
equatorial case. Nevertheless, it is clear that the essential features have been captured without the computationally 
intensive evaluation of the mode-sum out to 'infinity', or in reality, some radius very far from the horizon. 



X. CONCLUSIONS 



With the exception of Minkowski spacetime threaded by a cosmic string, explicit calculation of physically important 
renormalized expectation values in the axially symmetric geometry has proved elusive, though it is worth mentioning 
that Davies and Sahni [17] have calculated the vacuum polarization on the horizon of a Schwarzschild black hole 
threaded by a cosmic string for the restricted range (l/a) an integer, a result which we generalized to all a in [10]. 
Our calculation in this paper represents the first renormalization of the vacuum polarization on the exterior region 
of an axially symmetric black hole spacetime. Most importantly, we have renormalized without using properties 
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FIG. 5: Plot of {87vMy{{ip')ren ^ {^'}1T"''1 




\2/,-.2\oi^l 



FIG. 6: Plot of {8n MY {{^')ren - (r)*/™"') and (87rM)'^{<^-^)?,-;: 



of the special functions, (such as Addition Theorems for the Legendre functions, for example), but rather we have 
obtained useful summation formulae over the required mode functions that express geometrical singularities as non- 
convergent mode-sums. Moreover, we have appealed to the Hadamard singularity structure on a dimensionally reduced 
spacetime in order to guarantee the convergence of our renormalized expectation value. Our resultant expressions 
were rapidly convergent and relatively easy to compute numerically. This approach makes this method attractive in 
the renormalization process on the Kerr black hole, where we do not have the luxury of Addition Theorems. 

In addition, we have elucidated some points on how the direction in which we choose to point-split affects the order 
of our summation. We showed that a temporal splitting enforces an inner /, rn-sum and while for azimuthal splitting 
an inner n-sum is required. To the best of our knowledge, this subtlety has not previously been noted, perhaps 
because it has most often been convenient to split in the temporal direction in the spherically symmetric case; we 
became aware of it having, without due care in this direction, obtained different results for the same expression using 
azimuthal and temporal separation. Understanding of such issues is of key importance to extend our results to the 
Kerr space-time, since provisional calculations suggest that it is more convenient in that case to split in the azimuthal 
direction, as we have done in this paper. 
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Appendix 



In this appendix, we illustrate how to find the Hadamard divergence of the I, m mode-sum of the Green's function. 
We first note that since we are only interested in the divergent part, the order of summation is irrelevant since different 
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orders result in a finite discrepancy. Separating out the temporal dependence by writing 



G{x,x')^^ £ e-«(---')G„(x,x'), (11.1) 



Aw 

n— — oo 



and substituting into the wave equation, we find 



9,9 , d 1 d . ^ d 

smO 



\or or 



sine 09 89 a^sm^9d<l)^ 

r^ — 2Mr/ asmt^ 

Multiplying across by (1 — 2M/r)^^, this equation is a 3-dimensional Laplace equation with potential 

(V2 - y)G„(x, x') - -53"'/'<5(x - x') (11.3) 

where V^ is the invariant Laplacian on the 3-nietric 

Asl = dr^ + (r^ - 2Mr)d6|2 + (r^ - 2Mr)a^ siv? 9d(t)^ (11.4) 

where the potential is given by 

T) /<" T* 

An important result from Hadamard [18] is that the singularity structure of the Green's function is independent of 
the potential and depends only on the principal part of the wave operator. In 3 dimensions, the Hadamard form is 

^"("'"')^(2^^^ + ^("'"') (^^•') 

where 2ct is the square of the geodesic distance on the 3-mctric (11.4) and t/(x,x'), Vt^(x, x') are regular biscalars. 
One can expand U and W in terms of a [19], where to the required order we have 

t/(x,x') = l + 0{<j) 
T4^(x,x') =0(1). (11.7) 

For azimuthal separation, we have 

2cr = (r^ - 2Mr)a^ sin^ 9A4)^ + 0(A(/))4 (11.8) 

so that the Green's function at point separated in this way is 

^"(^' ^' ^^) - (.^-2M.)ksin^A0 + ^(1) 

^ 0(1). (11.9) 



M{f]^-iy/^asm9A(f> 



We also have the equivalent mode-sum expression for G„, obtained by comparing Eq. (11.1) with Eq. (2.11). Equating 
this mode-sum with expression (11.9) gives 

l—O m— — l \ I / 

In addition, wc have derived elsewhere [10] the following result for these Legendre Functions 

[2(l-cosa7)]i/2 -aZ^^Z.^^ r(A - |m|/a + 1)^^ (cosf^ji^^ (cos^j 



27ra 7o [2(1 - cos 9 cos 6*' + sin 6* sin^' coshu)] V2 
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where 

cosq;7 = cos 6* cos 6*' + sin sin 0' cosaA0 (11.12) 

and 

sm(* -7r/a) sm(* + 7r/a) m Ti\ 

l'a(u,^>) ^ — ^- — ^— — ^- — ^— . (11.13) 

cosn(u/Q;) — cos(W — TT/a) cosn(u/Q;) — cos(W + Tr/a) 

This formula relates a geometrical Hadamard singularity to a non-convergent mode-sum over the Legendre functions, 
valid for 1/2 < a < 1. The integral term here is regular everywhere (apart from the poles where there is a curvature 
singularity) and vanishes for a — I. For small (/) separations, this becomes 



1 



-E E -^""^ m^!1/"t!! ^A-'""°(cos^)- + o(i). (11.14) 

a -^-^ ^-^ 1 (A — m /q; + 1) 



aA^sin6' a. ^-^ ^-^_ r(A - |to|/q; + 1) 

Dividing across by M(j]^ — 1)^/^ and substituting into Eq.(ll.lO) yields 

Therefore, we have shown that the subtraction term here captures the divergence of the 3D Green's function G„ as 
a mode-sum. 

In a completely analagous way, one can associate the divergence of the /, m mode-sum in the subtraction terms of 
Eq. (6.4) with the divergence of 3D Green's function with a particular potential. The key point here is that this 3D 
divergence structure is captured entirely by the n — term since all n 7^ terms differ only by a potential and so do 
not affect the singularity [18]. Thus, for a given n, the singular part of the /, m sum of Eq. (6.4) is captured by the 
n = term 

y- y- »mA0(2A + l) ^(^ + 1"^'/" + ^^ ^-''"'/"(cosg)^ ^ fll 16) 

{r'^^J ^^ ^T{\-\m\/a + l)^ ^ ^ M2(r;-K 1)2 (2A + 1) ^ ' 

Comparing this to the subtraction term of Eq. (11.15) we see that we can guarantee the regularity of the Z, m sum 
for a given n by choosing ^ according to: 

iv + 1)2 
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